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1  Introduction 


The  Rayleigh  criterion  states  that  the  limit  of  angular  resolution  of  an  antenna  aperture  of 
size  D  is  6  ~  in  which  6  is  a  multiple  of  ® .  depending  on  the  geometry  of  aperture, 
and  A  is  the  wavelength  corresponding  to  the  transmitted  frequency.  The  Multiple  Signal 
Classification  (MUSIC)  (1)  algorithm,  using  the  eigenvector  decomposition  method,  is  a 
super-resolution  algorithm  widely  used  to  locate  closely  spaced  multiple  emitters  (targets) 
with  high  resolution  (better  than  the  Rayleigh  criterion).  The  algorithm  is  often  employed 
with  a  uniform  linear  array  of  sensors  of  total  aperture  length  D.  and  it  relies  on  the  phase 
progression  between  sensors  of  the  returned  signal.  In  contrast,  a  single  antenna,  also  of 
aperture  D  and  having  a  beam  width  9,  which  is  scanned  in  angle,  generates  no  phase  shifts 
within  the  main  lobe  because  of  the  target  direction;  consequently,  the  MUSIC  algorithm 
has  not  been  used  by  others  to  super-resolve  target  positions  within  the  main  beam  of  a 
single  narrow  beam  antenna.  In  this  report,  a  novel  technique  of  radar  array  processing  that 
employs  a  scanning  antenna  to  improve  target  detection  in  angular  domain  is  described  and 
assessed. 

The  algorithm,  termed  the  Scan-MUSIC  (Multiple  Signal  Classification)  or  SMUSIC,  was 
developed  by  the  U.S.  Army  Research  Laboratory  (ARL)  to  improve  the  spatial  resolution 
of  two  closely  spaced  point  targets  with  a  step-scanning  radar  antenna.  These  targets  were 
physically  located  within  the  beam  width  of  the  antenna.  However,  there  were  some  lim¬ 
itations  of  the  algorithm  for  two,  single  point  scatterer  targets  based  on  the  experimental 
and  computer  simulation  results  shown  in  references  (2)  through  (5).  The  algorithm  re¬ 
solved  the  targets  when  they  exhibited  constructive  interference;  it  could  not  resolve  the 
targets  when  they  exhibited  destructive  interference.  In  this  report,  the  application  of  the 
algorithm  was  extended  to  a  realistic  complex  scatterer-type  target  problem.  The  complex 
targets  used  were  located  within  the  beam  width  of  a  JLa-band  radar  antenna  with  frequency 
diversification,  and  the  targets  had  different  orientations  and  offsets. 

This  report  also  includes  a  graphical  interface  tool  to  run  and  analyze  the  algorithm  for  re¬ 
alistic  complex  scatterer-type  targets.  This  tool,  written  in  MATLAB,  enables  us  to  analyze 
the  perfomance  of  the  Scan-MUSIC  algorithm  in  some  practical  scenarios  and  is  a  design 
tool  to  assist  in  design  and  development  of  future  Army  radar  systems. 

The  report  is  divided  into  four  sections.  Section  2  describes  experimental  data  collection. 
Section  3  presents  the  raw  data,  the  SMUSIC  results,  and  the  tool  to  analyze  experimental 
data.  Section  4  draws  a  conclusion.  In  addition,  all  MATLAB®*  codes  of  the  tool  are 
included  in  appendix  A. 

*  MATLAB®  is  a  registered  trademark  of  The  MathWorks. 
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2  Experimental  Data  Collection 

Extensive  field  experimental  data  have  been  accumulated  over  the  years  with  a  linearly 
spaced  sensor  array.  However,  this  has  not  been  done  for  a  scanning  radar  antenna.  The 
objective  of  the  experiment  described  in  this  report  is  to  use  a  super-resolution  algorithm, 
SMUSIC,  for  the  analysis  of  the  measured  held  data  of  the  reflections  of  two  real,  complex 
scatterer-type  targets  via  a  scanning  antenna.  Subsection  2.1  of  this  section  describes  the 
experimental  setup.  Subsection  2.2  describes  the  data  collection.  Subsection  2.3  describes  a 
procedure  how  to  correct  the  phase  information  from  the  experimental  data  set. 

2.1  Experimental  Setup 

The  instrumentation  radar  system  operates  at  Ka  band  (35  GHz)  and  was  placed  on  the 
roof  top  of  a  truck.  A  circular  antenna  was  mounted  on  a  pedestal  and  had  an  aperture  of 
5  in.  with  conical  narrow  azimuth  beam  width.  The  radar  was  computer  controlled  and  had 
a  high  data  rate  recording  capability.  Figure  1  shows  the  AA-band  radar  mounted  on  the 
pedestal. 


Figure  1:  ATa-band  radar  and  pedestal. 

The  experimental  data  for  this  SMUSIC  analysis  were  collected  at  Aberdeen  Proving  Ground 
(APG),  MD.  The  test  range  was  a  hat,  grassy,  open  held.  Two  targets  (armored  personnel 
carrier  [APC]  M113  and  M60  tank)  were  placed  on  the  held  at  a  range  of  approximately 
600  m  from  the  radar.  These  two  targets  were  separated  about  1/3  of  the  antenna  beam 
width.  Table  1  shows  four  test  setups  for  this  experiment.  Figure  2  shows  the  hrst  setup 
(Run  7)  of  the  targets  that  were  placed  at  about  600  m  down-range  position  facing  toward 
the  radar.  The  range  was  measured  to  the  approximate  centroid  of  the  targets.  When  we 
look  at  the  hgure,  the  left  target  is  an  Ml  13  and  the  right  one  is  an  M60. 

Run  8,  shown  in  hgure  3,  kept  the  tanks  at  the  same  range  as  Run  7  but  rotated  the  M60 
tank  at  an  angle  approximately  40°. 
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Figure  2:  Run  7,  both  tanks  at  the  same  down-range  distance  measured  from  center  of  the 
radar  to  the  approximate  target  centroids. 


Figure  3:  Run  8,  M60  tank  rotated  approximately  40°. 
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Table  1:  Experimental  test  matrix  for  SMUSIC. 


Run  No. 

Figure 

Range  @  600  m 

Target  angle 

7 

2 

Same 

Head  on 

8 

3 

Same 

M60  40°  Offset 

9 

4 

5  m  closer 

Head  on 

10 

5 

5  m  closer 

M60  40°  Offset 

Run  9,  shown  in  figure  4,  left  the  position  of  the  M113  unchanged  but  moved  the  M60  tank 
to  the  radar  about  5  m  from  the  approximate  centroid  of  the  Ml  13  to  approximate  centroid 
of  the  M60  tank. 


Figure  4:  Run  9,  M60  tanks  moved  toward  to  the  radar  about  5  m  measured  from  the  two 
target  centroids  with  no  rotation. 

Run  10,  shown  in  figure  5,  maintained  the  same  target  location  as  Run  9  and  rotated  the 
M60  tank  at  the  same  angle  as  in  Run  8. 

The  antenna  amplitude  pattern  was  obtained  with  an  aluminum  trihedral  radar  cross  section 
(RCS)  =  19  dBsm.  The  reflector  was  placed  100  nr  down  range  on  the  radar  line-of-sight 
axis.  The  average  of  the  antenna  magnitude  data  over  the  frequency  band  was  used  for  the 
calculation  with  SMUSIC  algorithm. 
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Figure  5:  Run  10,  same  as  in  the  third  target  setup  with  the  M60  tank  rotated  approximately 
40°. 

2.2  Data  Collection 

The  signals  were  collected  as  the  antenna  step-scanned  horizontally  across  the  target  vehicles. 
At  each  step  angle,  the  radar  transmitted  a  coherent  stepped  frequency  waveform  which 
had  320  frequency  steps  from  32.4  GHz  to  35.6  GHz  with  10-MHz  increments  and  at  four 
different  polarizations  (vertical/vertical  [VV],  horizontal/horizontal  [HH],  vertical/horizontal 
[VH]  and  horizontal/vertical  [HV]).  The  in-phase  and  quadrature-phase  (I  &  Q)  channel  data 
were  stored  for  five  different  range  gates.  Figure  6  shows  the  layout  of  the  raw  data  stored 
as  integer  voltages  (from  analog-to-digital  converter)  in  variable  “out”.  In  this  figure,  RGi 
stands  for  range  gate  i,  where  i  =  1,  ...  5;  Az  angles  stand  for  azimuthal  step  angles. 
Frequency  (320)  is  the  index  of  stepped  frequency. 

The  following  is  the  MATLAB  code  used  to  format  data  as  shown  in  figure  6: 
data=reshape (out , [4  5  2  320  402]); 

in  which  "data"  is  raw  data  format;  "out"  is  a  variable  of  an  array 
of  5,145,600  elements  that  will  be  reformatted  (with  a  reshape  MATLAB 
function) in  the  format  of  4  receiver  channel  (VI,  VQ,  HI,  and  HQ),  5  range 
gates, 2  polarizations  (V  or  H) ,  320  step  frequencies,  and  402  step  angles. 


3  Results 

This  section  describes  the  raw  data  format,  the  description  of  SMUSIC  algorithm,  and  the 
performance  of  the  algorithm  and  a  tool  to  analyze  the  experimental  data. 
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RG1 


HH  (l,Q) 


Figure  6:  Raw  data  format. 
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3.1  Raw  Data 


The  raw  data  were  collected  and  stored  in  five  range  gates  to  verify  that  the  targets  were 
correctly  placed  in  range  gate  3.  Since  we  know  that  the  targets  were  placed  at  range  gate 
3,  we  use  range  gate  3  data  and  VV  polarization  for  data  processing.  The  motivation  of 
this  experiment  was  to  demonstrate  the  use  of  frequency  diversification  to  resolve  targets  in 
cross  range,  regardless  of  interferences.  Since  we  have  320  frequencies,  we  will  average  all 
the  frequencies  to  obtain  the  input  data  vector,  a?,  for  the  SMUSIC  algorithm.  The  raw  data 
format  is  shown  in  figure  7  for  range  gate  3. 


RG3 


Figure  7:  Range  gate  3  raw  data. 

Figures  8  to  11  are  magnitude  data  that  were  extracted  from  range  gate  3  for  each  run 
corresponding  to  targets  set  up  as  shown  in  figures  2  through  5  in  subsection  2.1.  We 
use  only  the  magnitude  data  in  this  case  because  the  phase  information  from  data  does 
not  add  much  value  in  this  calculation  because  of  the  noisy  and  “uncorrected”  phase  from 
the  experimental  data.  From  those  figures,  we  see  that  the  returned  signals  do  not  show 
how  many  target  centroids  there  are  in  the  responses  even  though  the  angular  separation 
between  two  targets  was  about  1/3  of  the  antenna  beam  width.  Our  goal  is  to  use  the 
SMUSIC  algorithm  to  resolve  target  centroids  from  those  data. 


3.2  Scan-MUSIC  Algorithm 

Super-resolution  algorithms,  such  as  the  MUSIC  algorithm,  have  a  special  structure  that 
uses  the  phase  augmentation  in  elements  of  a  linearly  spaced  sensor  array.  However,  for 
a  scanning  antenna,  the  phase  is  almost  constant  within  the  main  lobe  of  the  transmitted 
antenna  beam  width,  while  the  amplitude  of  the  response  varies.  Therefore,  the  returned 
signal  does  not  have  the  same  structure  as  for  an  array  of  uniformly  spaced  sensors. 


The  algorithms  that  use  the  signal  subspace  method  use  the  noise  eigenvector  subspace  from 

,  in  which  x  is 


xx 


H 


a  data  correlation  matrix.  The  correlation  matrix  is  defined  as  R  =  E 
the  input  data  vector  and  E  is  the  ensemble  average  from  the  snapshots.  Since  we  have  only 
one  snapshot  (one  complete  scan)  for  each  experimental  run,  the  SMUSIC  algorithm  fails 
to  resolve  closely  spaced  targets  within  the  main  lobe  of  the  antenna  beam  width  because 
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Figure  8:  Experiment  data,  Run  7,  range  gate  3,  polarization  =  VV. 
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Figure  9:  Experiment  data,  Run  8,  range  gate  3,  polarization  =  VV. 
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Figure  10:  Experiment  data,  Run  9,  range  gate  3,  polarization  =  VV. 
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Figure  11:  Experiment  data,  Run  10,  range  gate  3,  polarization  =  VV. 
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of  rank  deficiency  in  the  correlation  matrix.  To  remedy  this  and  to  “decorrelate”  coherent 
signals,  Evans  et  ah  (9)  introduced  the  spatial  smoothing  method,  and  its  proof  was  shown 
in  Shan  et  ah  (10)  and  Friedlander  et  ah  (11)-  This  method  was  originally  derived  for 
the  linearly  spaced  sensor  array.  We  adapted  the  spatial  smoothing  method  in  the  SMUSIC 
algorithm  for  a  scanning  radar  antenna,  called  the  subvector  averaging  method  because  the 
returned  signal  was  received  in  a  vector  as  the  antenna  stepped  through  discrete  angles  in 
the  angular  region  of  interest.  This  is  a  crucial  step  in  producing  a  correlation  matrix  of 
sufficient  rank  as  discussed  more  fully  in  reference  (2).  Following  is  the  description  of  the 
subvector  averaging  method. 

Given  input  data  x  as  described  in  section  3.1,  let  us  partition  the  input  data  vector,  x  = 
[xi,X2,  ■  ■  ■  ,Xn]t ,  with  N  step  scans  into  overlapping  subvectors  of  size  Q,  with  elements 


{1, . . . ,  Q}  forming  the  first  subvector,  with  elements  {2, . . 
subvector,  and  so  on,  as  follows: 

. ,  Q  +  1}  forming  the  second 

X{1}  =  [Xl  ,X2,-..,Xq}T, 

X{2}  =  [x2,  X3,  .  .  .  ,  Xq+\]T , 

* 

x{p}  =  [xp,xp+i, . . .  ,xn]t, 

(1) 

in  which  P  =  N  —  Q  +  1  is  the  total  number  of  subvectors,  P  <  Q. 

Next,  let  the  instantaneous  sub  correlation  matrix  R'1''1  be  x^x^H.  Finally,  define  the 
estimated  correlation  matrix  Rxx  to  be  the  average  of  P  subcorrelation  matrices  as 

ifaz  (2) 

^  p=  1  ^  p=  1 

The  matrix  Rxx  is  of  the  order  0  x  0 .  If  there  are  two  target  centroids,  there  will  be  two  pre¬ 
dominant  eigenvalues  with  the  orders  of  magnitude  larger  than  other  eigenvalues  from  Rxx- 
Let  {A,  E}  be  the  eigendecomposition  of  Rxx-,  where  the  eigenvalues  A  =  jAi,  A2, . . . ,  Aq| 
are  arranged  in  decreasing  order  of  magnitude,  with  corresponding  normalized  eigenvectors 
E  =  {ei,  e2, . .  • ,  6q}.  Let  us  partition  these  into  signal  and  noise  subspace  components  with 
Eg  =  {ei,  e2, . . . ,  €!m}  having  the  dimension  QxM  and  En  =  {&M+ i,  &m+ 2,  •  •  • ,  having 
the  dimension  Q  x  (Q  —  M ),  so  that  E  =  {Es,  EN}  and,  similarly,  Xs  =  {Ai,  A2, . . . ,  Am} 

and  Ajv  =  |Am+i,  Am+2,  •  •  • ,  Aq},  so  that  A  =  {As,  Xn}-  These  eigenvector  subspaces  are 
orthogonal  to  each  other.  For  this  report,  we  know  there  are  two  targets,  so  we  can  split  the 
signal  subspace  with  dimension  Q  x  2  and  the  noise  subspace  with  dimension  0  x  (Q-  2). 
In  a  general  case,  one  could  use  the  minimum  description  length  (MDL)  technique  (12,13) 
to  estimate  the  number  of  targets  before  splitting  the  signal  and  noise  subspaces. 
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Let  us  introduce  the  a(/3)  vector  to  be  the  effective  array  manifold  vector  for  all  allowable 
angles  (3  in  the  true  signal  space.  In  this  case,  a(/3)  is  the  two-way  magnitude  pattern  of 
the  /C,-band  radar.  Note  that  the  number  of  elements  of  the  manifold  vector  must  be  the 
same  number  of  elements  in  each  column  of  the  noise  eigenvector  subspace.  Since  the  signal 
subspace  contains  information  about  the  angles  of  arrival  from  each  plane  wave,  the  array 
manifold  vectors  from  those  angles  are  also  orthogonal  to  the  vectors  in  the  noise  subspace. 
The  magnitude  of  the  product  of  the  estimated  manifold  vector  and  the  noise  subspace  is 

therefore  a  small  number,  i.e.,  \\aH  {f3)E^\\2  =  e  or 

a((3)HENE*a((3)  =  e,  (3) 


in  which  e  is  a  small  number. 

Equation  3  is  a  null  SMUSIC  spectrum.  If  we  normalize  this  equation,  it  becomes 


a(0)HENE*a(P) 

a{(3)Ha((3) 


<  1. 


(4) 


The  inverse  of  equation  4,  i.e.,  the  magnitude  product  between  a  manifold  vector  from  all 
possible  angles  and  the  noise  subspace  matrix,  called  the  SMUSIC  pseudo-spectrum,  is 


P MU  {(3) 


oQ3)HaC3) 

a(f})»ENEHNa(li)' 


for  all  (3. 


(5) 


The  quantity  of  equation  5  is  large  for  [3  so  that  a((3) s  are  in  the  true  signal  space.  We 

see  that  the  peaks  from  this  pseudo-spectrum  are  more  intuitively  appealing.  Pmu((3)  is  the 
profile  of  a  one-dimensional  pseudo-spectrum.  The  locations  of  the  peaks  of  this  spectrum 
are  the  estimated  locations  of  the  target  centriods.  Equation  5  only  estimates  the  locations, 
not  their  strength.  There  is  a  separate  task  to  determine  the  strength  of  the  estimated 
signals  by  application  of  the  least  square  solution  technique  (1).  This  report  does  not  cover 
this  task. 

3.3  Numerical  Implementation  of  the  SMUSIC  Algorithm 

Equation  5  has  the  same  form  as  an  ordinary  MUSIC  calculation.  The  SMUSIC  algorithm 
performs  the  search  for  all  possible  angles  [3  in  the  searching  angular  region  to  locate  the 
peaks  of  the  SMUSIC  pseudo-spectrum.  This  searching  method,  which  goes  through  each 
computational  loop  as  shown  in  figure  12  for  each  searching  angle,  consumes  a  lot  of  com¬ 
putational  time.  A  new  procedure  or  method  has  been  developed,  which  takes  advantage  of 
the  optimized  MATLAB  software,  for  a  fast  computation  of  the  SMUSIC  spectrum  as  shown 
in  figure  13.  The  technique  involves  the  following: 
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1.  Pre-calculate  a(/3)  at  each  angle  for  all  the  searched  angles. 

2.  Store  them  in  a  hie  or  a  computer  memory. 

3.  Create  a  searched  array  manifold  matrix,  <A{f3)  in  which  the  columns  of  the  matrix  are 
the  array  manifold  vector  of  the  searched  angles. 


Simulated  data 


.  Partition  input  vector 
.  Form  submatrices 
.Average  submatrices 


Using  SVD 


Using  SVA 


Compute 

Determine  signal 

eigenvalues  and 

- ► 

subspace  and  noise 

eigenvectors 

_ 

subspace 

\  4 

Determining  £^and  E 


N 


A  fR)  - 

k 

Locate 

peaks  of  SMUSIC  spectrum 

MP)f   «Op)«(P) 

MU  IP  / 

\ 

h"«(  p) 

2  «ff(p  )£*£"<?(  p) 

\ 

Figure  12:  Ordinary  SMUSIC  process  for  experimental  data  and  simulated  data. 
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The  SMUSIC  pseudo-spectrum  is  shown  below  by  the  use  of  the  property  of  the  matrix 
product  and  the  diagonal  matrix  manipulation  to  develop  this  fast  technique: 


fast  MU 


(«  = 


Aqg(-4(/j)H-4(/?)) 

diag(A(p)HENEHNAm 


for  all  p. 


(6) 


From  equation  6,  we  see  that  the  numerator  is  the  vector  of  all  the  diagonal  elements  of 
the  dot  product  of  the  searched  array  manifold  matrix  and  its  conjugate  transpose.  Since 
we  know  «A(/3),  the  only  thing  that  we  need  in  the  denominator  is  to  compute  the  noise 
eigenvector  subspace  EN. 


Figure  13:  Fast  SMUSIC  process  for  experimental  data  and  simulated  data. 

In  summary,  from  those  figures,  SVA  stands  for  subvector  averaging  method,  and  SVD  stands 
for  singular  value  decomposition.  By  comparing  the  ordinary  MUSIC  as  shown  in  figure  12 
and  its  fast  implementation  in  figure  13,  we  see  that  the  equation  in  the  last  block  of  the 
figure  12  is  replaced  with  equation  6  for  the  fast  implementation  as  shown  in  figure  13.  We 
used  the  fast  implementation  for  the  analysis  in  the  simulation.  The  fast  implementation  is 
about  100  times  faster  than  the  ordinary  implementation. 
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3.4  Tool  to  Analyze  Experimental  Data 

In  this  section,  we  describe  an  interface  program.  The  interface  program  was  written  in 
the  MATLAB  language,  version  6  or  R12.  The  MATLAB  code  for  this  program  is  listed  in 
appendix  A.  To  run  this  program,  first  you  need  to  invoke  the  MATLAB  command  window. 
Then,  assumming  that  the  interface  program  is  in  the  current  directory,  type  smusicinter 
Ver3  at  the  MATLAB  prompt  sign  in  the  MATLAB  command  window.  The  SMUSIC 
interface  window  will  appear  on  the  screen.  Figure  14  shows  the  interface  program  window 
which  consists  of  the  plotting  window  and  the  input  parameter  section.  The  input  parameter 
section  includes  three  parts:  input  hie,  processing,  and  plotting  options. 


Figure  14:  Interface  program  window. 
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3.4.1  Input  File  Selection 

This  selection  allows  users  to  select  different  experimental  data.  The  data  hies  used  are 

Run7,  Run8,  Run9,  and  RunlO,  as  described  section  2. 

3.4.2  Processing 

This  section  allows  users  to  select 

•  Range  Gate  -  five  range  gates  (1  to  5). 

•  Polarization  -  VV,  VH,  HV,  HH. 

•  Comp. /Mag.  -  compute  SMUSIC  pseudo-spectrum  using  either  complex  or  magnitude 
data. 

•  FREQ  -  select  a  particular  frequency  or  the  average  of  all  frequencies  for  the  SMUSIC 
calculation. 

•  SVA  -  Number  of  subvector  averages,  i.e.,  the  number  of  times  that  we  perform  subvec¬ 
tor  averaging  to  estimate  the  data  correlation  matrix.  For  example,  SVA  =  10  means 
that  we  perform  10  partitions  of  the  input  data,  form  these  subcorrelation  matrices, 
and  then  calculate  the  average  of  those  subcorrelation  matrices. 

•  Run  SMUSIC  -  the  execution  button  starts  the  SMUSIC  calculation. 

3.4.3  Plotting  and  Other  Options 

This  section  allows  users  to  select 

•  Plot  SMUSIC  -  plots  SMUSIC  pseudo-spectrum  for  the  simulation  after  the  SMUSIC 
calculation  is  done. 

•  Plot  TGTSETUP  -  plots  experimental  target  setup,  depending  on  the  input  hie 
selection. 

•  Plot  M.  Ave.  Data  -  plots  average  magnitude  of  experimental  data.  This  calculation 
is  the  average  of  all  frequency  steps  within  the  radar  beam  width. 

•  Plot  P.  Ave.  Data  -  plots  average  phase  of  experimental  data.  This  calculation  is  the 
average  of  all  frequency  steps  within  the  radar  beam  width. 

•  M.  @  freq  1  ...  -  plots  magnitude  raw  data  at  a  paticular  frequency  which  is  used  for 
the  SMUSIC  calculation. 

•  P.  @  freq  1  ...  -  plots  phase  raw  data  at  a  particular  frequency  which  is  used  for  the 
SMUSIC  calculation. 

•  Plot  M.  Ave.  Ant.  -  plots  average  magnitude  of  the  measured  antenna  pattern.  This 
calculation  is  the  average  of  all  frequency  steps  within  the  radar  beam  width. 
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•  Plot  P.  Ave.  Ant.  -  plots  average  phase  of  the  measured  antenna  pattern.  This 
calculation  is  the  average  of  all  frequency  steps  within  the  radar  beam  width. 

•  M.  @  freq  1  ...  -  plots  magnitude  antenna  at  a  particular  frequency  which  is  used  for 
the  SMUSIC  calculation. 

•  P.  @  freq  1  ...  -  plots  phase  at  a  particular  frequency  which  is  used  for  the  SMUSIC 
calculation. 

•  Info  -  This  option  shows  where  the  MATLAB  source  code  of  the  interface  program  is 
located. 

•  ZOOM  -  is  an  option  to  allow  users  to  span  the  figure  at  a  specific  area.  To  do  this, 
first  click  the  ZOOM  button,  and  then  move  the  mouse  cursor  in  the  desired  area  of 
the  figure  and  click  the  left  mouse  button.  If  you  want  to  zoom  back  to  the  normal 
size  figure,  click  the  right  mouse  button. 

•  REFRESH  -  clears  all  input  parameters,  figure,  and  screen  command  window. 

•  QUIT  -  exits  the  interface  program. 

3.5  SMUSIC  and  Experimental  Data  Results 

Figures  15  through  18  show  the  windows  of  the  interface  program  and  the  results  of  the 
SMUSIC  algorithm  with  selected  parameters  for  the  simulation.  We  selected  experimental 
data  at  range  gate  3,  VV  polarization,  average  magnitude  data.  These  figures  show  the 
performance  of  the  SMUSIC  algorithm  and  the  experimental  data  as  a  function  of  azimuth 
angles.  In  each  figure,  the  solid  curve  is  the  SMUSIC  result  and  the  dash-dot-dash  curve  is  the 
experimental  data.  The  shaded  area  is  the  — 3-dB  beam  width  of  the  antenna.  The  “peaks” 
of  the  solid  curve  are  the  estimation  of  the  target  centroid  locations.  We  see  clearly  that  we 
achieve  our  goal  to  resolve  the  targets  that  are  well  within  the  antenna  beam  width.  In  fact, 
the  separation  of  the  two  peaks  is  about  1/3  of  antenna  beam  width  for  four  experimental 
runs. 

In  summary,  the  SMUSIC  algorithm  resolved  the  targets  that  were  placed  within  the  — 3-dB 
beam  width  of  the  iQ-band  radar  antenna.  Although  we  collected  both  real  (in-phase)  and 
imaginary  (quadrature)  experimental  data,  we  used  only  the  magnitude  data  that  were  the 
square  root  of  the  sum  of  square  of  the  real  and  the  imaginary  to  perform  the  SMUSIC 
calculation.  The  combination  of  the  real  and  imaginary  data  formed  complex  experimental 
data. 


4  Conclusions 

The  report  has  described  the  novel  extension  of  the  MUSIC  algorithm  to  the  Scan-MUSIC 
(SMUSIC)  for  complex  scatterer-type  targets.  The  extension  showed  clearly  the  processing 
scheme  of  the  SMUSIC  algorithm  which  used  a  single  rotatable  antenna  without  increasing 
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Figure  15:  SMUSIC  result  and  magnitude  experimental  data  for  Run  7,  range  3,  Polarization 
=  VV,  SVA  =  10. 


Figure  16:  SMUSIC  result  and  magnitude  experimental  data  for  Run  8,  range  3,  Polarization 
=  VV,  SVA  =  10. 
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Figure  17:  SMUSIC  result  and  magnitude  experimental  data  for  Run  9,  range  3,  Polarization 
=  VV,  SVA  =  10. 


Figure  18:  SMUSIC  result  and  magnitude  experimental  data  for  Run  10,  range  3,  Polariza¬ 
tion  =  VV,  SVA  =  10. 
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an  antenna  aperture.  We  used  the  SMUSIC  algorithm  with  magnitude  experimental  data  to 
resolve  two  closely  spaced  targets  (an  M60  tank  and  M113  APC)  and  to  estimate  the  locations 
of  their  centroids  in  cross  range  with  frequency  diversification,  regardless  of  orientation 
and  offset  of  the  two  targets.  We  show  that  we  can  improve  the  angular  resolution  of  the 
targets  with  SNR  about  10  to  12  dBsm  by  about  a  factor  of  3.  We  also  described  the 
graphical  user  interface  program  which  allows  us  to  analyze  the  performance  of  the  SMUSIC 
algorithm.  This  report  also  described  a  fast  implementation  of  the  SMUSIC  algorithm,  which 
shows  great  promise  in  being  implemented  in  a  real-time  radar  system.  The  success  of  this 
implementation  will  enhance  future  U.S.  Army  radar  systems. 
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A  MATLAB  source  code  of  the  Interface  Program 

7,  Filename:  smuexpinter_Ver3 .m  (SMUSIC  experimental  interface  file) 

7,  Author :  Canh  Ly 
7.  Date :  18  August  2003 
7. 

7.  SMUSICEXPINTER_Ver2  plots  experimental  data  and  SMUSIC  process  using  GUI. 

7.  Example : 

7.  Input  File:  smusic_run07 

7.  Processing:  Range  Gate:  Range  Gate  3 

7.  Polarization:  VV 

7.  Comp.  /Mag.  :  Complex 

7.  One/Two  Way :  Two  way 

7.  SVA  :  10 

7.  Run  SMUSIC 

7. 

7.  Ver  2:  18  Aug  2003,  Deleted  One/Two  Way  Control  Button 
/  Added  FREQ  to  select  indidual  frequency  or  the 

7.  average  of  all  frequencies 

7.  Ver  3:  9  February  2004:  Delete  Phase  (Uncorrected  and  Corrected)  button  only  used  uncorrected 

nsteps  =  l;value=l;  7.  Setup  default  number  of  steps  and  value  for  frequency  selection 
nfreqs  =  1; 
phase_value=l ; 

runno=7;range_gate=3;pol=’VV’ ; 

figtag  =  findobj(0,  ’tag1,  ’SMUSIC  -  Experimental  Complex  Target  (Ver.  3)’); 
if  isempty (figtag) 

fig  =  f igure ( ’Name ’ ,  ’SMUSIC  -  Experimental  Complex  Target  (Ver.  3) ’, ’Units’ , ’pixels’ , 
’Position’,  [310  250  700  495],... 

’NumberTitle’ ,  ’off’,  ’resize’,  ’on’,... 

’tag’,  ’SMUSIC  -  Experimental  Complex  Target  (Ver.  3)’); 

end 

set(gcf,  ’Def aultAxesPosition’ ,  [.075  .15  .6  .75]); 

figcolor  =  get(gcf,  ’color’); 

hndtexthead. Style  =  ’text’; 
hndtexthead. Units  =  ’normal’; 
hndtexthead. backgroundcolor  =  figcolor; 
hndtexthead. foregroundcolor  =  ’blue’; 

hndt ext subhead. Style  =  ’text’; 

hndt ext subhead. Units  =  ’normal’; 

hndt ext subhead. backgroundcolor  =  figcolor; 

hndt ext subhead. foregroundcolor  =  ’black’; 

7.  Input  File  Text  position 
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uicontrol  (gcf ,  hnd.texth.ead,... 

’position1,  [.755  .90  .15  .0385],... 

’ f ontsize ’  ,  10 , ’ f ontweight ’ , ’bold’ , . . . 

’string’,  ’Input  File’); 

7,  Processing  Text  position 
uicontrol (gcf ,  hndtexthead, . . . 

’position’,  [.755  .80  .15  .0385],... 

’ f ontsize ’ , 10 , ’ f ontweight ’ , ’bold’ , . . . 

’string’,  ’Processing’); 

uicontrol (gcf ,  hndtext subhead, ’Horiz’ , ’center’ , 
’position’,  [.72  .76  .115  .035],... 

’string’,  ’Range  Gate’); 

uicontrol (gcf ,  hndtext subhead, ’Horiz’ , ’center’ , 
’position’,  [.873  .76  .115  .035],... 

’string’,  ’Polarization’); 

uicontrol (gcf ,  hndtext subhead, ’Horiz’ , ’center’ , 
’position’,  [.72  .68  .115  .035],... 

’string’,  ’Comp. /Mag. ’) ; 

7,  uicontrol (gcf ,  hndtext subhead, ’Horiz ’,’ center 
7.  ’position’,  [.873  .68  .115  .035],... 

7.  ’string’,  ’One/Two  Way’); 

uicontrol (gcf ,  hndtext subhead, ’Horiz’ , ’center’ , 
’position’,  [.873  .68  .115  .035],... 

’ string’ ,  ’FREQ ’ ) ; 

7.  uicontrol  (gcf  ,  hndtextsubhead, ’Horiz ’,’ center 
7.  ’position’,  [.72  .595  .10  .035],... 

7.  ’string’,  ’Phase’); 

7.  uicontrol (gcf ,  hndtextsubhead, ’Horiz ’,’ center 
7.  ’position’,  [.873  .595  .10  .035],... 

7.  ’string’,  ’SVA’); 

uicontrol (gcf ,  hndtextsubhead, ’Horiz’ , ’center’ , 
’position’,  [.80  .595  .10  .035],... 

’ string’ ,  ’ SVA’ ) ; 


7.  Editable  buttons 
hndedit . Style  =  ’edit’; 
hndedit .Units  =  ’normal’; 
hndedit .Font size  =  10; 
hndedit . f ontweight  =  ’bold’; 

7.  Define  popup  button 
hndpopup . Style  =  ’popupmenu’ ; 
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hndpopup. Units  =  'normal’; 

7.  Define  push,  button 
hndpush. Style  =  ’push’; 
hndpush. Units  =  ’normal’; 

inputfilestr=[’Run  7’; 

’Run  8’ ; 

’Run  9’ ; 

’Run  10’] ; 

inputfile=  uicontrol (gcf ,  hndpopup,  ’position’,  [.76  .85  .134  .05],... 
’string’ ,  inputf ilestr , ’backgroundcolor ’ , ’white’ , ’enable’ , ’on’ , . . . 
’Interruptible’ , ’on’ , ’Horiz’ , ’center’ ,  ... 

’foregroundcolor’ , ’black’) ; 

set (inputf ile , ’ callback’ , ’ setf ile ’ ) ; 

rangestr=[’ Range  Gate  1 1  Range  Gate  2  I  Range  Gate  3  I  Range  Gate  4  I  Range  Gate  5’] 
rangeinput  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.71  .72  .134  .05],... 
’string’ ,  rangestr, ’backgroundcolor’ , ’white’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’black’) ; 

set (rangeinput ,  ’callBack’ , ’ setrangegate ’ ) ; 

polstr=[’V  VlH  HIV  H I H  V’] ; 

polinput  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.87  .72  .115  .05],... 

’ string ’ ,  polstr , ’ backgroundcolor ’ , ’ white ’ , ’ enable ’ , ’ on ’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’black’) ; 

set (polinput ,  ’ callBack’ , ’ setpol ’ ) ; 


cmstr= [ ’ Complex  I  Magnitude ’ ] ; 

cminput  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.71  .635  .134  .05],... 
’string’ ,  cmstr, ’backgroundcolor’ , ’white’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’black’) ; 

set (cminput ,  ’ callBack’ , ’ setprocess ’ ) ; 

7,  otstr=[’Two  way’]; 

7.  otinput=  uicontrol (gcf ,  hndpopup,  ’position’,  [.87  .635  .115  .05],... 
7.  ’string’,  otstr, ’backgroundcolor ’, ’white’ , ’enable’ , ’on’ ,.. . 

7.  ’  Interruptible ’,’ on’ ,  ... 

7.  ’foregroundcolor ’, ’black’ )  ; 

7.  set  (otinput ,  ’  callback’  ,  ’  setonetwo  ’  )  ; 

[data,ds]  =  extract21Apr03(runno,range_gate,pol) ; 
nfreqs=size (ds . IQ_data, 1) ; 
labelstrm= [] ; 
if  nfreqs  ==  1 

labelstrm=[’Freq.  1|’]; 
elseif  nfreqs>=2 
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temp  =  [’Freq.  1 | * ] ; 
for  i  =  2:nfreqs 

tempadd  =  [’Freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  freqs’]; 

sfreqdata  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.87  .635  .115  .05],.. 
’string’ ,  labelstrm, ’backgroundcolor ’ , ’white’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’black’) ; 
freqsel  =  get (geo , ’value ’) ; 

set (sfreqdata, ’callback’ ,  ’ selectfreqdata’ ) ; 


mrammrammrarammrafflmmmmmmiimm 


7,  19  Mar  2003:  Take  out  One  way  option  out 


lot str= [ ’ One  way  I  Two  way ’ ] ; 

70ot input =  uicontrol (gcf ,  hndpopup,  ’position’,  [.77  .46  .115  .05],... 
7, ’string’ ,  otstr, ’backgroundcolor’ , ’white’ , ’enable’ , ’on’ , . . . 

7«’ Interruptible ’, ’on’ ,  ... 

7. ’foregroundcolor ’  ,  ’black’)  ; 

7.set  (otinput ,  ’  callback’  ,  ’  setonetwo  ’  )  ; 


7.  otstr=[’Two  way’]; 

7.  otinput=  uicontrol (gcf ,  hndpopup,  ’position’,  [.77  .46  .115  .05],... 
7.  ’string’,  otstr, ’backgroundcolor ’, ’white’ , ’enable’ , ’on’ ,.. . 

/  ’ Interruptible ’,’ on’ ,  ... 

7.  ’foregroundcolor ’, ’black’ )  ; 

7.  set  (otinput ,  ’  callback’  ,  ’  setonetwo  ’  )  ; 


7. 

7.  phasestr=[’ Uncorrected  I  Corrected’]  ; 

7.  phaseinput  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.71  .555  .134  .05] 
7.  ’string’ ,  phasestr, ’backgroundcolor’ , ’white’ , ’enable’ , ’on’ , . . . 

7.  ’  Interruptible ’,’ on’ ,  ... 

7.  ’foregroundcolor ’, ’black’ )  ; 

7.  set  (phaseinput ,  ’  callBack’  ,  ’  setphase  ’  )  ; 

nniniiin  9  Feb  2004 


phase_value=l ; 
if  phase_value==l 

phastr=’ Uncorrected  phase’; 
elseif  phase_value==2 

phastr=’ Corrected  phase’; 

end 


7.  7.  Input  SVA 

7.  inputl  =  uicontrol (gcf ,  hndedit,  ’position’,  [.87  .56  .115  .045]); 
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7,  calll  =  [’sva  =  str2num(get  (geo ,  ’  ’string’  ’));’]  ; 

7,  set(inputl,  ’callBack’,  calll); 

inninniin  9  February  2004  mrararararamray.rarararay. 


7.  Input  SVA 

inputl  =  uicontrol (gef ,  hndedit ,  ’position’,  [.79  .56  .115  .045]); 
calll  =  [’sva  =  str2num(get (geo ,  ’’string’’));’]; 
set (inputl,  ’callBack’,  calll); 


7.  Run  SMUSIC 

calcstart  =  uicontrol (gef ,  hndpush,  ’position’,  [.79  .495  .115  .05],.. 
’string’ ,  ’Run  SMUSIC’ , ’backgroundcolor ’ , ’red’ , . . . 

’foregroundcolor’ , ’white’ , ’callback’ ,  ’ smusic_mf rf _inter_Ver2 ’ ) ; 

7.  Plot  SMUSIC  and  TGTSETUP 

plotsmu  =  uicontrol (gef ,  hndpush,  ’position’,  [.71  .43  .13  .05],... 
’string’,  ’Plot  SMUSIC’ , ’backgroundcolor’ , ’blue’ ,.. . 

’foregroundcolor’ , ’yellow’ , ’callback’ ,  ’plotsmuresult ’ ) ; 

plottgt  =  uicontrol (gef ,  hndpush,  ’position’,  [.855  .43  .13  .05],... 
’string’,  ’Plot  TGTSETUP’ , ’backgroundcolor’ , ’blue’ ,.. . 
’foregroundcolor’ , ’yellow’ , ’callback’ ,  ’plottgtsetup’ ) ; 

7.  Plot  data 

7.  First  row,  Left  button 

plotmdata  =  uicontrol (gef ,  hndpush,  ’position’,  [.71  .36  .13  .05],... 
’string’,  ’Plot  M.  Ave  Data’ , ’backgroundcolor’ , ’cyan’ ,.. . 
’foregroundcolor’ , ’black’ , ’callback’ ,  ’plotmaveexpdata’ ) ; 

7.  First  row,  right  button 

plotpdata  =  uicontrol (gef ,  hndpush,  ’position’,  [.855  .36  .13  .05],... 
’string’,  ’Plot  P.  Ave  Data’ , ’backgroundcolor’ , ’cyan’ ,.. . 
’foregroundcolor’ , ’black’ , ’callback’ ,  ’plotpaveexpdata’ ) ; 


7.  Second  row,  left  button 
7.  Initialized  string  for  magnitude  plot 


labelstrm= [] ; 
if  nfreqs  ==  1 

labelstrm= [’M.  @  freq.  1 | ’ ] ; 
elseif  nfreqs>=2 

temp  =  [’M.  @  freq.  1 | ’ ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’M.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  M.  data’]; 

plotmdatafreq  =  uicontrol (gef ,  hndpopup,  ’position’,  [.71  .30  .13  .05] 
’string’ ,  labelstrm, ’backgroundcolor’ , ’cyan’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 
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’ foregroundcolor ’ , ’black1) ; 
column  =  get (geo , ’value ’) ; 

set (plotmdatafreq, ’callback’ ,  ’plotmf reqdata’ ) ; 


7,  Second  row,  right  button 
7,  Initialized  string  for  phase  plot 
labelstrp  =  []  ; 
if  nfreqs  ==  1 

labelstrp=[’P.  @  freq.  1 1 ’ ]  ; 
elseif  nfreqs>=2 

temp  =  [’P.  @  freq.  1 | ’ ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’P.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrp  =  [temp  tempadd] ; 
temp  =  labelstrp; 

end 

end 

labelstrp= [labelstrp  ’Average  P.  data’]; 

plotpdatafreq  =  uicontrol (gef ,  hndpopup,  ’position’,  [.855  .30  .13  .05] 
’string’ ,  labelstrp, ’backgroundcolor ’ , ’cyan’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’black’) ; 
column  =  get (geo, ’value’) ; 

set (plotpdatafreq, ’ callback’ ,  ’plotpf reqdata’ ) ; 

7.  Plot  antenna 

7.  First  row,  left  button 

plotmaveant  =  uicontrol (gef ,  hndpush,  ’position’,  [.71  .23  .13  .05],... 
’string’,  ’Plot  M.  Ave  Ant .’, ’backgroundcolor ’, ’magenta’ ,.. . 
’foregroundcolor’ , ’white ’ , ’ callback’ ,  ’plotantpat ’ ) ; 

7.  First  row,  right  button 

plotpaveant  =  uicontrol (gef ,  hndpush,  ’position’,  [.855  .23  .13  .05],.. 
’string’,  ’Plot  P.  Ave  Ant .’, ’backgroundcolor ’, ’magenta’ ,.. . 
’foregroundcolor’ , ’white ’ , ’ callback’ ,  ’plotantpat ’ ) ; 

7.  Second  row,  left  button 
7.  Initialized  string  for  magnitude  plot 
labelstrm= [] ; 
if  nfreqs  ==  1 

labelstrm= [’M.  @  freq.  1 | ’ ] ; 
elseif  nfreqs>=2 

temp  =  [’M.  @  freq.  1 | * ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’M.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  M.  Ant.’]; 
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plotmantf req  =  uicontrol(gcf ,  hndpopup,  ’position’,  [.71  .17  .13  .05], 
’string’ ,  labelstrm, ’backgroundcolor ’ , ’magenta’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’white’) ; 
column  =  get (geo , ’value ’) ; 

set (plotmantf req, ’callback’ ,  ’plotmfreqant ’ ) ; 

7.  Second  row,  right  button 
7,  Initialized  string  for  magnitude  plot 
labelstrm= [] ; 
if  nfreqs  ==  1 

labelstrm=[’P.  @  freq.  1 1 ’ ] ; 
elseif  nfreqs>=2 

temp  =  [’P.  @  freq.  1 | ’ ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’P.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  P.  Ant.’]; 

plotpantfreq  =  uicontrol(gcf ,  hndpopup,  ’position’,  [.855  .17  .13  .05] 
’string’ ,  labelstrm, ’backgroundcolor’ , ’magenta’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’white’) ; 
column  =  get (geo , ’value ’) ; 

set (plotpantfreq, ’callback’ ,  ’plotpfreqant ’ ) ; 


7.  Push  bottons 

uicontrol (gef ,  hndpush,  ’position’,  [.775  .095  .10  .05],... 
’string’,  ’ INFO ’,’ callback’ ,  ’infofile’); 

uicontrol (gef ,  hndpush,  ’position’,  [.685  .025  .08  .05],... 
’string’,  ’ZOOM’ , ’callback’ ,  ’zoom’); 

uicontrol (gef ,  hndpush,  ’position’,  [.775  .025  .10  .05],... 
’string’,  ’REFRESH’ ,’ callback’ ,  ’ smuexpinterinit_Ver2 ’ ) ; 

uicontrol (gef ,  hndpush,  ’position’,  [.885  .025  .08  .05],... 
’string’,  ’QUIT’ , ’callback’ ,  ’close  all’); 
grid 
clc ; 
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7.  function  setfile 
7,  Filename:  setfile. m 
7,  Author  :  Canh  Ly 
7.  Date  :  16  March  2001 

file_value  =  get (geo, ’value’) ; 
if  f ile_value==l  7.  FFT  calculation 
titlestr=’Run  7’; 
runno=7 ; 

A=imread( ’ run7setup.jpg’ , ’ jpg’) ; 
elseif  f ile_value==2 
titlestr=’Run  8’; 
runno=8 ; 

A=imread( ’run8setup.jpg’ , ’jpg’) ; 
elseif  f ile_value==3 
titlestr=’Run  9’; 
runno=9 ; 

A=imread( ’run9setup.jpg’ , ’jpg’) ; 
elseif  f ile_value==4 
titlestr=’Run  10’; 
runno=10; 

A=imread( ’runlOsetup .jpg’ , ’ jpg’ ) 

end 

7.  function  setrangegate 
7.  Filename:  setrangegate. m 
7.  Author  :  Canh  Ly 
7.  Date  :  10  December  2001 

range_gate  =  get (geo, ’value’) ; 


7.  function  setpol 

7.  Filename:  setpol. m 

7.  Author  :  Canh  Ly 

7.  Date  :  10  December  2001 

pol_value  =  get (geo, ’value’) ; 

if  pol_value==l  '/,  FFT  calculation 
pol=’VV’ ; 

elseif  pol_value==2 
pol=’HH’ ; 

elseif  pol_value==3 
pol=’VH’ ; 

elseif  pol_value==4 
pol=’HV’ ; 

end 


7«function  setprocess 
7.  Filename:  setprocess. m 


7,  Author  :  Canh  Ly 

7,  Date  :  16  January  2001 

comp_mag  =  get (geo, ’value’) ; 

7,f unction  selectfreqdata(freqsel) 

7.  Filename:  plotmf reqdata.m 
7.  Author  :  Canh  Ly 
7,  Date  :  18  August  2003 
7. 

7.  Purpose:  Select  data  at  particular  frequency  for  SMUSIC  processing. 

7. 


freqsel  =  get (geo , ’value 1 ) ; 

7.  function  setphase 
7,  Filename:  setphase. m 
7.  Author  :  Canh  Ly 
7.  Date  :  14  August  2003 

phase_value  =  get (geo , ’value 1 ) ; 
if  phase_value==l 

phastr=’ Uncorrected  phase1; 
elseif  phase_value==2 

phastr=’ Corrected  phase’; 

end 


29 


7  Filename:  smusic_mfrf_inter_Ver2.m 
7,  Author :  Canh  Ly 
7.  Date :  8  September  2003 

7.  1  April  2003:  Modified  the  search  algorithm  as  a  fast  search  algorithm  using  the 
7.  property  of  the  diagonal  matrix  operation 

az= [42. 0:0. 2: 44. 8  45.0:0.01:55.0  55.02:0.2:57.0]; 

load  amvmatout46  7.  This  matrix  file  has  a  variable  Vout,  generated  from 
7.  genamvmat.m  25  April  2003 


rararamraramrararararararanrara 

7.  Open  antenna  file  and  calculate  phase  correction 

mmmmmmmmmmmmmmnmm 

[ant5,cal5]  =  antpatl8Apr03 (5 , ’VV’ , [1 :320] ) ;  7.  All  frequencies 


77777.7777.7.  start  from  here:  Old  procedure  of  correcting  phase777.777.77.7777.7.77 
7. 


7.  tempant=ant5.IQ_data_raw. ’ ; %  tempant:  320x791;  ant5.IQ_data_raw:  791x320 
7  ffttempant=f ft (tempant) ; 

7.  magf  fttempant=abs  (f  fttempant)  ; 

7.  phaf fttempant=angle (f fttempant)  ; 


7.  Calculate  the  phase  at  the  maximum  location  of  the  average  of  all  frequencies 
7.  [Y ,  I]  =max  (ant5 .  yantdB_raw)  ; 

7.  f  f  tantmax=f  ft  (ant5 .  IQ_data_raw(I ,  :))  ; 

7.  [maxant ,  Imaxant]  =max  (abs  (f  f  tantmax)  )  ; 

7.  phasef  f  tantmax=angle  (f  f  tantmax  (Imaxant))  ; 


7.  [magffttempantmax, indexff ttempantmax] =max(20*logl0(magffttempant) ) ; 

7.  anglef f ttempant=angle (f fttempant (indexff ttempantmax) ) ; 

7.  anglef fttempoff set=anglef fttempant -phasef f tantmax; 

7.  anglef fttempoff set=anglef f ttempant-anglef fttempant (1) ;  7.  Correct  all  angle  with 

7.  repect  to  the  first  location 

7.  A=anglef  fttempoff  set  ’  ;  7.  791x1  dimension 
7.  B=A(  :,  ones  (1 , 320)  )  ;  7.  791x320  dimension 
7.  C=B  .  ’  ; 

7.  tempD=phaf fttempant  (  :  ,  232  : 682)  -C  (  :  ,  232 :  682)  ; 

7.  D= [phaf fttempant (:, 1 : 231)  tempD  phaf fttempant (:, 683 : 791)] ; 

7.  E=D .  ’  ; 

7. 


7.  correctantmag=abs  (ant5 .  IQ_data_raw)  ; 

7.  correctantpha=E; 

7.  correctantpat=correctantmag.  *exp(j  .*E)  ; 
7. 


7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.77.7.7.7.7. 


End  of  old  procedure  of  correcting  phase 


7.7.7.77.7.7.7. 


7.7.7.77.7.777.7.  Start  from  here:  new  procedure  of  correcting  phase 
7.  19  August  2003 
7. 

ant5 . IQ_data_raw_un=ant5 . IQ_data_raw; 


77.7.777.7.7. 777.7.7.77. 
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if  phase_value==2 

phaseerror  7,  Calculate  phase  error  from  phaseerror.m 
phaerrormat=phaerror ( : , ones (1 , size (ant 5 . IQ_data_raw, 1) ) ) ; 
antphaerror=unwrap (angle (ant5 . IQ_data_raw . ’ ) ) -phaerrormat ; 
magIQ_data=abs (ant5 . IQ_data_raw) ; 
phaIQ_data=antphaerror . ’ ; 

ant5 . IQ_data_raw_cor=magIQ_data. *exp(j . *phaIQ_data) ; 
ant5 . IQ_data_raw  =  ant5 . IQ_data_raw_cor ; 
elseif  phase_value==l 

ant5 . IQ_data_raw  =  ant5 . IQ_data_raw_un; 

end 


mmm  nmmmmm 


End  of  new  procedure  of 


correcting  phase  TLVLVLTL 


7,  Initialized  string  for  antenna  magnitude  plot 
nfreqs=size (ant5 . IQ_data_raw,2) ; 
labelstrm= [] ; 
if  nfreqs  ==  1 

labelstrm=[’M.  @  freq.  1|’]; 
elseif  nfreqs>=2 

temp  =  [’M.  @  freq.  1 | * ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’M.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  M.  Ant.’]; 

plotmantfreq  =  uicontrol(gcf ,  hndpopup,  ’position’,  [.71  .17  .13  .05],... 
’string’ ,  labelstrm, ’backgroundcolor ’ , ’magenta’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 

’foregroundcolor’ , ’white’) ; 
column  =  get (geo , ’value ’) ; 

set (plotmantfreq, ’callback’ ,  ’plotmfreqant ’ ) ; 


7,  Initialized  string  for  antenna  phase  plot 
abelstrm=  [] ; 
if  nfreqs  ==  1 

labelstrm=[’P.  @  freq.  1 1 ’ ] ; 
elseif  nfreqs>=2 

temp  =  [’P.  @  freq.  1 | ’ ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’P.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  P.  Ant.’]; 

plotpantfreq  =  uicontrol(gcf ,  hndpopup,  ’position’,  [.855  .17  .13  .05],... 
’string’ ,  labelstrm, ’backgroundcolor’ , ’magenta’ , ’enable’ , ’on’ , . . . 

’ Interruptible ’,’ on’ ,  ... 
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’ foregroundcolor ’ , ’white1) ; 
column  =  get (geo , ’value ’) ; 

set (plotpantf req, ’callback’ ,  ’plotpfreqant ’ ) ; 


nmfflmrafflrarafflrannrafflfflffl 


7,  Open  data  file 


7,[data,ds]  =  extract (runno, range_gate, pol)  ; 
[data,ds]  =  extract21Apr03 (runno, range_gate, pol) ; 
nfreqs=size (ds . IQ_data, 1) ; 


7,  Added  on  18  Aug  2003  to  correct  phase  of  input  data 
ds . IQ_data_old=ds . IQ_data; 
if  phase_value==2 

phaseerror  7.  Calculate  phase  error  from  calibration  file 
pha_dataerrormat=phaerror ( : , ones (1 , size(ds . IQ_data, 2) ) ) ; 
phadsIQdata=unwrap (angle (ds . IQ_data) ) -pha_dataerrormat ; 
magdsIQdata=abs (ds . IQ_data) ; 

ds.IQ_data  =  magdsIQdata. *exp(j . *phadsIQdata) ; 
elseif  phase_value==l 

ds.IQ_data  =  ds . IQ_data_old; 

end 


7.  Initialized  string  for  data  magnitude  plot 

labelstrm  =  []  ; 

hndpopup . Style  =  ’popupmenu’ ; 

hndpopup. Units  =  ’normal’; 

if  nfreqs  ==  1 

labelstrm=[’M.  @  freq.  1 1 ’ ]  ; 
elseif  nfreqs>=2 

temp  =  [’M.  @  freq.  1 | * ] ; 
for  i  =  2: nfreqs 

tempadd  =  [’M.  @  freq.  ’  num2str(i)  ’I’]; 
labelstrm  =  [temp  tempadd] ; 
temp  =  labelstrm; 

end 

end 

labelstrm= [labelstrm  ’Average  M.  data’]; 

plotmdatafreq  =  uicontrol (gef ,  hndpopup,  ’position’ , [.71  .30  .13  .05],... 
’string’ ,  labelstrm, ’backgroundcolor ’ , ’cyan’ , ’enable’ , ’on’ , . . . 
’Interruptible’ , ’on’ ,  ... 

’foregroundcolor’ , ’black’) ; 
column  =  get (geo , ’value ’) ; 

set (plotmdatafreq, ’callback’ ,  ’plotmf reqdata’ ) ; 

7.  Initialized  string  for  data  phase  plot 
labelstrp  =  []  ; 
if  nfreqs  ==  1 

labelstrp=[’P.  @  freq.  1 | ’ ] ; 
elseif  nfreqs>=2 

temp  =  [’P.  @  freq.  1 | ’ ] ; 
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for  i  =  2:nfreqs 

tempadd  =  [ ’ P .  @  freq.  ’  num2str(i)  ’I’]; 
labelstrp  =  [temp  tempadd] ; 
temp  =  labelstrp; 

end 

end 

labelstrp= [labelstrp  ’Average  P.  data1]; 

plotpdatafreq  =  uicontrol (gcf ,  hndpopup,  ’position’,  [.855  .30  .13  .05], 
’string’ ,  labelstrp, ’backgroundcolor ’ , ’cyan’ , ’enable’ , ’on’ , . . . 
’Interruptible’ , ’on’ ,  ... 

’foregroundcolor’ , ’black’) ; 
column  =  get (geo, ’value’) ; 

set (plotpdatafreq, ’ callback’ ,  ’plotpf reqdata’ ) ; 

if  freqsel  ==  size(ds . IQ_data, 1)+1 

inputdata=[ds.outdatadB-min(ds.outdatadB)] . ’ ; 

else 

if  comp_mag==l 

inputdata=ds . IQ_data (freqsel , :) .  ’  ; 
else 

inputdata=abs (ds . IQ_data(f reqsel ,:)).’; 

end 

end 

Rxhat=estrxhatfs (inputdata, sva)  ; 

[eigvalO , eigvecO] =eigvv(Rxhat)  ; 

[signalO , noiseO] =snspaces (eigvalO , eigvecO , 2) ; 

[nrownoisesum,ncolsnoisesum] =size(noiseO) ; 

esumO=noiseO; 

sumsuba_sum=0 . 0 ; 

ncolA=size (Vout , 2) ; 

suba_sum  =  Vout(l :nrownoisesum, : ) ; 

num  =  diag(suba_sum’ *suba_sum) ; 

dem  =  diag(suba_sum’ * (esumO*esumO ’ ) *suba_sum) ; 

pmusics0=real(10*logl0(num. /dem) ) ; 

plot  (az  ,pmusics0,  ’b’  )  ;hold 
7,  set(gca, ’Xlim’ ,  [40  58]) 
if  (max(pmusicsO)>=  max (inputdata) ) 
axis ([35  60  0  max(pmusics0)+5] ) ; 
else 

axis ([35  60  0  max (inputdata) +5] ) ; 

end 

plottitlestr=[titlestr, ’ ,  Range  gate  =  ’ ,num2str (range_gate) ,  ... 

’,  Polarization  =  ’,pol,’,  5"  in  ant.  @  35  GHz’]; 
title(plottitlestr) ; 
xlabel(’ Azimuth  angle  (deg)’) 
grid  on 
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function  [ant, cal]  =  antpatl8Apr03(runno,pol,freqs) 

7,  Filename:  antpatl8Apr03 .m 
7,  Author :  Canh  Ly 
7.  Date:  18  April  2003 
7. 

7.  Purpose:  The  program  extracts  data  from  smusic_run5 . sl6  files  from  the  field 
7.  test  on  25  October  2001. Then  it  plots  the  antenna  pattern  in  which 

7.  the  corner  reflector  was  placed  at  100  m  (one  range  gate  at  11)  with 

7.  the  RCS  =  19.0  dBsm,  or  79.4  m~2. 

7. 

7.  Syntax : 

7.  [ant, cal]  =  antpatl8Apr03(runno,pol,freqs) ; 

7.  where : 

7.  ant  :  data  structure,  output  antenna  pattern. 

7. 

7.  Note : 

7.  Data  for  the  antenna  pattern  is  a  multidimensional  array.  The  dimension  is 
7.  shown  below : 

7.  -  first  dimension  is  receiver  channel,  4  (VI,  VQ,  HI  and  HQ) 

7.  -  second  dimension  is  range  gate,  1  (r=l) 

l  -  third  dimension  is  polarization,  2  (p=  V  or  H) 

7.  -  fourth  dimension  is  number  of  step  frequencies,  320 

7.  -  fifth  dimension  is  azimuthal  angles,  1002 

7.  The  total  size  of  the  data  file  is  =  4*1*2*320*1002  =  2565120 

7.  Example : 

7.  1.  Run  antenna  file  number  5,  VV  polarization  at  freqs  190  to  290  or  34.4  to  35.4  GHz 

7.  [ant5 , cal5]  =  antpatl8Apr03 (5, ’VV’ ,  [190:290]  )  ; 

7.  [ant6 , cal6]  =  antpatl8Apr03 (6,  ’VV’  ,  [190:290]  )  ; 

7.  2.  Run  antenna  file  number  5,  VV  polarization  at  freqs  1  to  320  or  33.2  to  35.4  GHz 

7.  [ant5 ,  cal5]  =  antpatl8Apr03(5 ,  ’  VV’  ,  [1 :  320]  )  ; 

7.  [ant6,cal6]  =  antpatl8Apr03(6 , ’ VV’ , [1 : 320] ) ; 

7.  3.  Run  antenna  file  number  5,  VV  polarization  at  freqs  100  to  320  or  34.2  to  35.4  GHz 

7.  [ant5 , cal5]  =  antpatl8Apr03(5,  ’VV’  ,  [100:320]  )  ; 

7.  [ant6 , cal6]  =  antpatl8Apr03 (6, ’VV’ ,  [100:320]  )  ; 

7. 

if  runno  ==  5 

sfilename  =  ’D:\cly\SMUSIC_MFRF\data\smusic_run5.sl6’ ; 
load  D : \cly\SMUSIC_MFRF\data\smusic_run5 . tag ; 

tagdata=smusic_run5 ;  7.  Used  to  extract  the  second  column  for  the  azimuthal  angles, 
elseif  runno==6 

sfilename  =  ’D:\cly\SMUSIC_MFRF\data\smusic_run6.sl6’; 
load  D : \cly\SMUSIC_MFRF\data\smusic_run6 . tag ; 

tagdata=smusic_run6 ;  7.  Used  to  extract  the  second  column  for  the  azimuthal  angles. 

end 

7.  Open  calibration  file 
[cal]  =  calibration(5) ; 

titlestr=’ Antenna  Pattern’; 

if  pol==’VV’  I  pol==’vv’  I  pol==’Vv’  I  pol==  ’vV’ 
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receiver_pol  = 

i; 

pol_selection= 

i; 

elseif  pol==,VHI  I 

pol==,vh’ 

1  pol==,VhI 

1  pol==  ’  vH 

receiver_pol  = 

3; 

pol_selection= 

i; 

elseif  pol==,HHI  I 

O 
1— 1 

II 

II 

1  pol==,HhI 

1  pol==  >hH 

receiver_pol  = 

3; 

pol_selection= 

2; 

elseif  pol==’HV’  I 

pol==’hv’ 

1  pol==’Hv’ 

1  pol==  >hV 

receiver_pol  = 

i; 

pol_selection= 

2; 

end 

titlestring= [titlestr  ’ ,  polarization  =  ’  pol] ; 

7,  Read  experimental  data  file  and  store  the  data  in  a  multidimensional  array, 
fid  =  f open(sf ilename , ’rb’ , ’b’) ; 
out  =  fread(fid, ’intl6’) ; 

data=reshape  (out ,  [4  1  2  320  1002]);  7,  See  the  above  comments 

7.  Extract  offset  data  at  the  first  azimuth  angle,  the  first  angle  of  402  angles. 
7.  Use  3.2  GHz  bandwidth,  corresponding  from  frequency  1  to  frequency  320. 
off set=data( : , : , : ,freqs, 1) ; 

tfs  =  length (freqs)  ;  7«  Total  number  of  frequency  steps 

7.  Calcualte  average  offset  for  VI  channel 
off set_VVI=of f set (1,1,1,:); 

off set_HVI=of f set (1,1,2,:);  7.  Added  on  7  April  2003 
7.  For  all  frequencies 

avg_of f set_VVI=sum(of f set_VVI ( : ) )/tf s ; 
avg_of f set_HVI=sum(of f set_HVI ( : ) )/tf s ; 
ant . avg_of f set_VVI  =  avg_of f set_VVI ; 

ant . avg_of f set_HVI  =  avg_of f set_HVI ;  7.  Added  on  7  April  2003 

7.  Calculate  average  offset  for  VQ  channel 
of f set_VVQ=of f set (2,1,1,:); 

off set_HVQ=of f set (2,1,2,:);  7.  Added  on  7  April  2003 
avg_of f set_VVQ=sum(of f set_VVQ( : ) )/tf s ; 
avg_of f set_HVQ=sum(of f set_HVQ( : ) )/tf s ; 
ant . avg_of f set_VVQ  =  avg_of f set_VVQ ; 

ant . avg_of f set_HVQ  =  avg_of f set_HVQ ;  7.  Added  on  7  April  2003 

7.  Calculate  average  offset  for  HI  channel 
of f set_VHI=of f set (3,1,1,:); 

off  set_HHI=of  f  set  (3,1,2,:);  7.  Added  on  7  April  2003 
avg_of f set_VHI=sum(of f set_VHI ( : ) )/tf s ; 
avg_of f set_HHI=sum(of f set_HHI ( : ) )/tf s ; 
ant . avg_of f set_VHI  =  avg_of f set_VHI ; 

ant . avg_of f set_HHI  =  avg_of f set_HHI ;  7.  Added  on  7  April  2003 
7.  Calculate  average  offset  for  HQ  channel 
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off set_VHQ=of f set (4,1,1,:); 

off  set_HHQ=of  f  set  (4,1,2,  :  )  ;  7.  Added  on  7  April  2003 
avg_of f set_VHQ=sum(of f set_VHQ( : ) )/tf s ; 
avg_of f set_HHQ=sum(of f set_HHQ( : ) )/tf s ; 
ant . avg_of f set_VHQ  =  avg_of f set_VHQ ; 

ant .  avg_of  f  set_HHQ  =  avg_of  f  set_HHQ ;  °i  Added  on  7  April  2003 
ant . avg_of f set_VV  =  [avg_of f set_VVI  avg_offset_VVQ] ;  7,  Added  on  7  April  2003 

ant . avg_of f set_VH  =  [avg_of f set_VHI  avg_offset_VHQ] ;  7«  Added  on  7  April  2003 

ant . avg_of f set_HV  =  [avg_of f set_HVI  avg_offset_HVQ] ;  %  Added  on  7  April  2003 

ant . avg_of f set_HH  =  [avg_of f set_HHI  avg_offset_HHQ] ;  7«  Added  on  7  April  2003 

7. 

7.  Take  out  the  second  repeated  data  value  in  the  signature  file. 

7.  First  append  one  arbitrary  number  at  the  beginning  of  the  file. 

7.  Then  find  the  difference  of  all  data,  if  the  same,  the  output  is  zero. 

7,  Detect  those  zeros,  and  use  only  non-zero  data. 

7. 

theta  =  tagdata(2 : size (tagdata, 1) ,2) ; 
nzeroindex=f ind( (diff ( [1 ;theta] ) ) ~=0) ; 
ant .theta_raw=theta(nzeroindex) ; 
theta  =  theta (nzeroindex) ; 

l  Extract  data 

7.  Extract  VI  data  for  320  frequencies  and  401  azimuthal  angles  at  range  gate  3 
7.  Use  the  full  bandwidth 

VIdata=reshape(data(receiver_pol , 1 ,pol_selection,freqs , 2 : 1002) ,tf s , 1001) ; 

7.  VIdata  after  taking  out  repeated  data  due  to  the  pedestal  jitter  because 
7.  of  small  anlge  step  size. 

VIdata=VIdata( : .nzeroindex) ; 

VQdata=reshape(data(receiver_pol+l , 1 ,pol_selection, freqs , 2 : 1002) , tf s , 1001) ; 

7.  VQdata  after  taking  out  repeated  data  due  to  the  pedestal  jitter  because 
7.  of  small  anlge  step  size. 

VQdata=VQdata( : .nzeroindex) ; 

7.  Since  there  were  two  receivers,  V  and  H  each  has  I  and  Q  (  or  VI ,  VQ,  HI  and  HQ), 
7,  there  were  two  independent  transmitters  V  and  H.  Therefore,  we  can  use  VV  (I  &  Q) 
7.  and  VH  (I  &  Q)  to  calculate  the  offsets.  Other  channels  such  as  HH  (I  &  Q)  and  VH 
7.  I  &  Q)  are  almost  identical  to  others. 

if  pol==’VV’  I  pol==’vv’  I  po^^Vv’  I  pol==  ’  vV1 
I_data=(VIdata-ant . avg_of f set_VVI) ; 

Q_data=(VQdata-ant . avg_of f set_VVQ) ; 
elseif  pol==,VHI  I  pol==’vhV  I  pol==,Vh’  I  pol==  ’  vH’ 

I_data=(VIdata-ant . avg_of f set_VHI) ; 

Q_data=(VQdata-ant . avg_of f set_VHQ) ; 
elseif  pol==’HH’  I  pol==’hh>  I  pol==>Hh’  I  pol==  >hH’ 

I_data=(VIdata-ant . avg_of f set_VHI) ; 

Q_data=(VQdata-ant . avg_of f set_VHQ) ; 
elseif  pol==’HV’  I  pol==’hv>  I  pol==’Hv’  I  pol==  >hV’ 

I_data=(VIdata-ant . avg_of f set_VVI) ; 

Q_data=(VQdata-ant . avg_of f set_VVQ) ; 

end 
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fflXffl 

7,  Calculate  the  average  of  the  magnitude  antenna  pattern  for  selected  band  of  frequency. 

xmm 

7,  For  all  frequencies 
outdata=sum(I_data. ~2+Q_data. ~2)/tfs; 
ant . I_data_raw=I_data’ ; 
ant . Q_data_raw=Q_data’ ; 

ant . IQ_data_raw=ant . I_data_raw+ j . *ant . Q_data_raw ; 

IQ_data=I_data+j . *Q_data; 
ant . yant_raw=outdata. ’ ; 

ant . yantdB_raw=(10*logl0(outdata) ) . ’ ;  7.  Use  10*logl0  because  outdata  is  the 

7.  average  of  the  sum  of  square 

7.  1 .  Compute  the  magnitude  response  of  the  antenna  response  at  each  frequency 

7.  2.  Determine  the  maximum  angle  based  on  1 .  at  all  frequencies  (101) 

7.  3.  Compute  the  average  peak  from  2. 

7.  4.  Compute  the  error  of  the  peak  of  antenna  response  at  each  frequency  with 

7.  respect  to  the  average  peak 

7.  5.  Calculate  the  standard  deviation  of  the  error  of  the  peak  of  all  antenna 
7.  responses. 

magantmatrix=abs (ant . IQ_data_raw) ; 

[maxmat ,maxmatindex] =max(magantmatrix) ; 
aveofpeak  =  sum(ant . theta_raw (maxmat index) ) /tf s ; 
errorpeak=ant . theta_r aw (maxmat index) -aveofpeak; 
stdpeak  =  std(errorpeak) ; 

7.  1 .  Locate  the  maximum  antenna  response 

7.  2.  Normalize  azimuthal  scan  angle  to  obtain  "symmetrical"  antenna 
[maxmag,maxmagindex] =max(max(abs (ant . IQ_data_raw) ) ) ; 

7of igure (1) ;plot (ant . theta_raw- (aveofpeak-stdpeak) , abs(ant . IQ_data_raw( : ,maxmagindex) ) ) ; 

outdatamag=abs (ant . IQ_data_raw( : ,maxmagindex) ) ; 

outdatacomplex=ant . IQ_data_raw( : ,maxmagindex) ; 

theta_f req81=  ant . theta_raw- (aveofpeak-stdpeak) -0 . 0844 ; 

[azscans81 ,magscans81 , comp81] =ant smooth (theta_f req81 , outdatamag,outdatacomplex, 15) ; 
skippts=3 ; 

ant.azscan81  =  azscans81 (1 : skippts : length(azscans81) ) ; 

ant.comp81  =  comp81 (1 : skippts : length(comp81) ) ; 

ant .magscan81=  (magscans81 (1 : skippts : length(magscans81) ) ) ; 

7.-min(magscans81  (1 :  skippts  :  length(magscans81)  )  )  ; 

[Y,I]  =  max(outdata) ; 

yant=(outdata-Y) . ’ ; 

yantdB= (10*logl0 (outdata) ) . ’ ; 

theta  =  theta-theta(I) ; 

ant. theta  =  theta; 

ant . yant=yant ; 

ant . yantdB=yantdB ; 

yant_test=(ant . yantdB_raw-min(ant ,yantdB_raw) ) ; 
y ant  _t  e  stmax=y ant  _t  e  s t -max ( y ant  _ t  e  s t ) ; 

7.  f  igure ; plot  (ant .  theta, yant_testmax,  ’b’  ) 
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if  runno==5 

thetaindex=f ind(ant . y ant dB -max (ant . yantdB)<=-2 . 92  &  ant .yantdB-max(ant .yantdB) >=-3 . 07)  ; 

angletwoindex=ant . theta (thetaindex) ; 

ant .beamwidth  =  angletwoindex(2)-angletwoindex(l) ; 

end 

lessmlOindex=f ind(theta<=-10 . 0) ; 

gtmlOindex=f ind(theta>-10 . 0)  ; 

thetapatch=- (theta(lessmlOindex)+l) ; 

thetaout  =  [theta (gtmlOindex) ;  flipud(thetapatch)] ; 

thetaindex=f ind(thetaout>=-10  &  thetaout<=10) ; 

xantout  =  thetaout (thetaindex) ; 

ant.thetap  =  xantout; 

yantpatch=yant (lessmlOindex) ; 

yanttemp  =  [yant (gtmlOindex) ;  flipud(yantpatch)] ; 
yantout  =  yanttemp (thetaindex) ; 

ant.yantp  =  yantout; 
yantdBpatch=yantdB (lessmlOindex) ; 

yantdBtemp  =  [yantdB (gtmlOindex) ;  flipud(yantdBpatch)] ; 
yantdBout  =  yantdBtemp (thetaindex) ; 

ant . yantdBp  =  yantdBout ; 

I_datapatch=I_dat a (lessmlOindex) ; 

I_datatemp=  [I_data(gtmlOindex) ;  flipud(I_datapatch)] ; 
ant . I_data=I_datatemp (thetaindex) ; 

Q_datapatch=Q_dat a (lessmlOindex) ; 

Q_datatemp=  [Q_data(gtmlOindex) ;  flipud(Q_datapatch)] ; 
ant . Q_data=Q_datatemp (thetaindex) ; 

IQ_datapatch=IQ_data (lessmlOindex) ; 

IQ_datatemp=  [IQ_data (gtmlOindex) ;  flipud(IQ_datapatch)] ; 
ant . IQ_data=IQ_datatemp (thetaindex) ; 

scanindex=f ind(ant . thetap>=-8 . 5  &  ant . thetap<=8 . 5) ; 
azscantemp=ant . thetap (scanindex) ; 
magscantemp=ant . yantdBp(scanindex) ; 

[azscans ,magscans , mags cans] =ant smooth (azscantemp,magscantemp,magscantemp, 3) ; 
skippts=2 ; 

ant.azscan  =  azscans (1 : skippts : length(azscans) ) ; 

ant .magscan=  (magscans (1 : skippts : length (mags cans) ) ) -min (mags cans (1 : skippts : length (mags cans) ) ) ; 
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function  [data,ds]  =  extract21Apr03(runno,range_gate,pol) 

7,  Filename:  extract21Apr03 .m 

7,  Author :  Canh  Ly 

7.  Date:  21  April  2003 

7. 


Purpose:  The  program  extracts  data  from  *.sl6  files  from  the  field  test  on 
25  October  2001. 

Then  it  plots  the  magnitude  response  at  the  ’ range_gate’  and  ’pol’ 
input  selection. 

Syntax : 

[data,ds]  =  extract21Apr03(runno,range_gate,pol) ; 
where : 

tagdata:  input  tagfile  indicating  how  many  azimuthal  angles  in  each  run. 
data:  output  data  is  a  multidimensional  array.  The  dimension  is  shown  below: 

1.  For  data  file  (not  antenna  pattern) 

-  first  dimension  is  receiver  channel,  4  (c=  I  or  Q) ,  VI,  VQ ,  HI  and  HQ 

-  second  dimension  is  range  gate,  5  (  r=  1 ,  2,  3,  4,  5) 

-  third  dimension  is  polarization,  2  (p=  V  or  H) 

-  fourth  dimension  is  number  of  step  frequencies,  320 

-  fifth  dimension  is  azimuthal  angles,  402 

The  total  size  of  the  data  file  is  =  4*5*2*320*402  =  5145600 

2.  For  antenna  pattern  file 

-  first  dimension  is  receiver  channel,  4  (c=  I  or  Q) ,  VI,  VQ ,  HI  and  HQ 

-  second  dimension  is  range  gate,  1  (r=l) 

-  third  dimension  is  polarization,  2  (p=  V  or  H) 

-  fourth  dimension  is  number  of  step  frequencies,  320 

-  fifth  dimension  is  azimuthal  angles,  1002 

The  total  size  of  the  data  file  is  =  4*1*2*320*1002  =  2565120 
ds  :  data  structure  for  the  average  of  the  offset  values  as  the  background 
runno :  run  number  from  data  signatures 

[7]  Run  number  7  for  SMUSIC_run7 

[8]  Run  number  8  for  SMUSIC_run8 

[9]  Run  number  9  for  SMUSIC_run9 

[10]  Run  number  10  for  SMUSIC_runl0 

range_gate:  range  gate  selection  (1  to  5) 

pol  :  polarization  (>W’,  >HH’,  ’VH> ,  or  >HV’) 

Example : 

[data,ds]  =  extract21Apr03(7,3, ’VV’) ; 


t  Modified:  8  April  2003:  Added  selection  of  frequency  from  190  to  290  because 
7.  of  flat  phase  response. 

7. 


if  runno  ==  7 

sf ilename  =  >D : \CLY\SMUSIC_MFRF\data\smusic_run7 . sl6 ’ ; 
load  D : \CLY\SMUSIC_MFRF\data\smusic_run7 . tag; 
tagdata=smusic_run7 ; 
titlestr=’Run  7’; 
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elseif  runno  ==  8 

sf ilename  =  ’D : \CLY\SMUSIC_MFRF\data\smusic_run8 . sl6 ’ ; 
load  D : \CLY\SMUSIC_MFRF\data\smusic_run8 . tag; 
tagdata=smusic_run8 ; 
titlestr=’Run  8’; 
elseif  runno  ==  9 

sf ilename  =  ’D : \CLY\SMUSIC_MFRF\data\smusic_run9 . sl6 ’ ; 
load  D : \CLY\SMUSIC_MFRF\data\smusic_run9 . tag; 
tagdata=smusic_run9 ; 
titlestr=’Run  9’; 
elseif  runno  ==  10 

sf ilename  =  ’D : \CLY\SMUSIC_MFRF\data\smusic_runlO . sl6 ’ ; 
load  D : \CLY\SMUSIC_MFRF\data\smusic_runlO . tag; 
tagdata=smusic_runlO ; 
titlestr=’Run  10’; 

end 


if  pol==’VV’  I  pol==’vv’  I  pol==’Vv’  I  pol==  ’vV’ 
receiver_pol  =  1; 
pol_selection=  1; 

elseif  pol==’VH’  I  pol==’vh’  I  pol==’Vh’  I  pol==  ’vH’ 
receiver_pol  =  3; 
pol_selection=  1; 

elseif  pol==’HH’  I  pol==’hh’  I  pol==’Hh’  I  pol==  ’hH’ 
receiver_pol  =  3; 
pol_selection=  2; 

elseif  pol==’HV’  I  pol==’hv’  I  pol==’Hv’  I  pol==  ’hV’ 
receiver_pol  =  1; 
pol_selection=  2; 

end 


titlestring= [titlestr  ’,  at  range  gate  ’  num2str (range_gate)  and  polarization  =  ’  poll; 

fid  =  fopen(sf ilename , ’rb’ , ’b’) ; 
out  =  fread(fid, ’intl6’) ; 

data=reshape  (out ,  [4  5  2  320  402]);  '/,  See  the  above  comments 

7,  Extract  offset  data  at  the  first  azimuth  angle,  the  first  angle  of  402  angles. 

7.  Only  use  2.2  GHz  bandwidth,  corresponding  to  from  frequency  100  to  frequency  320  (221  freqs) . 
7,  of  f  set=data(  100  : 320 , 1)  ; 

7.  Only  use  1.0  GHz  bandwidth,  corresponding  to  from  frequency  190  to  frequency  290  (101  freqs). 
7.  from  34.4  GHz  to  35.4  GHz  (8  April  2003) 

7.  of  f  set=data(  190  : 290 , 1)  ; 

off  set=data(  1 :  320 , 1)  ;  7.  For  all  frequencies 

7.  Modified:  31  Mar  2003 
7.  For  VI  channel 

7.  Compute  the  average  offset  of  VI  channel: 

7.  1.  Obtain  the  offset  of  all  range  gates  (1  to  5)  for  all  frequencies  (221  from 
7.  freq  100  to  320)  .  The  offset_VVI  has  a  5  x  221  dimension 
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7,  2.  Calculate  the  average  offset  of  VI  by  averaging  all  frequencies  (size(offset_VVI,2)) 

7,  all  columns)  first  and  then  for  all  range  gates  (size(off set_VVI , 1) )  (all  rows). 

7.  offset_VVI  =  reshape  (of  f  set  (1 1 ,:),  [5  221]); 

7.  8  April  2003:  Choose  frequencies  from  190  to  290 
7.  offset_VVI  =  reshape  (of  f  set  (1 1 ,:),  [5  101]); 

offset_VVI  =  reshape (of f set (1 1 ,:), [5  320]); 

ds . avg_of f set_VVI  =  sum (sum (of f set_VVI ,2) / size (of f set_VVI ,2) )/ (size (of f set_VVI , 1) ) ; 

7.  For  average  offset  of  VQ  channel 

7.  offset_VVQ  =  reshape  (of  f  set  (2 ,  :  ,  1 ,  :  )  ,  [5  221]);  For  frequencies  from  100  to  320 
7.  8  April  2003:  Choose  frequencies  from  190  to  290 
7.  offset_VVQ  =  reshape  (of  f  set  (2 ,  :  ,  1 ,  :  )  ,  [5  101]); 
offset_VVQ  =  reshape (of f set (2 , : , 1 , : ) , [5  320]); 

ds . avg_of f set_VVQ  =  sum (sum (of f set_VVQ ,2) / size (of f set_VVQ ,2) ) / (size (of f set_VVQ , 1) ) ; 

7.  For  HI  channel 

7.  offset_VHI  =  reshape  (of  f  set  (3,  :  ,  1 ,  :  )  ,  [5  221]);  '/,  For  frequencies  from  100  to  320 
7.  8  April  2003:  Choose  frequencies  from  190  to  290 
7.offset_VHI  =  reshape  (of  f  set  (3,  :  ,  1 ,  :  )  ,  [5  101]); 
offset_VHI  =  reshape (of f set (3, : , 1 , : ) , [5  320]); 

ds . avg_of f set_VHI  =  sum (sum (of f set_VHI ,2) / size (of f set_VHI ,2) )/ (size (of f set_VHI , 1) ) ; 

7.  For  HQ  channel 

7.  offset_VHQ  =  reshape  (of  f  set  (4,  :  ,  1 ,  :  )  ,  [5  221]);  For  frequencies  from  100  to  320 
7.  8  April  2003:  Choose  frequencies  from  190  to  290 
7.  offset_VHQ  =  reshape  (of  f  set  (4,  :  ,  1 ,  :  )  ,  [5  101]); 
offset_VHQ  =  reshape (of f set (4, : , 1 , : ) , [5  320]); 

ds . avg_of f set_VHQ  =  sum (sum (of f set_VHQ ,2) / size (of f set_VHQ ,2) )/ (size (of f set_VHQ , 1) ) ; 

7. 

7,  Take  the  repeated  data  values  out  in  the  signature  file. 

7.  First  append  one  arbitrary  number  at  the  beginning  of  the  file. 

7.  Then  find  the  difference  of  all  data,  if  the  same,  the  output  is  zero. 

7.  Detect  those  zeros,  and  use  only  non  zero  data. 

7. 

theta  =  tagdata(2 : size (tagdata, 1) ,2) ; 
ds.theta_raw  =  theta; 

nzeroindex=f  ind(  (diff  (  [1  ;theta]  )  )  ~=0)  ;  7.  Pick  out  no  repeated  data 
theta  =  theta (nzeroindex) ; 
l  Extract  data 

7.  Extract  VI  data  for  320  frequencies  and  401  azimuthal  angles  at  range  gate  3 
7.  Only  use  2.2  GHz  bandwidth,  corresponding  from  frequency  100  to  frequency  320  (221  freqs) . 
7.  VIdata=reshape (data(receiver_pol ,range_gate ,pol_selection, 100 : 320,2 : 402) ,221,401) ; 
VIdata=reshape(data(receiver_pol ,range_gate ,pol_selection, 1 : 320,2: 402) ,320,401) ; 

7.  VIdata  after  taking  out  repeated  data  due  to  the  pedestal  jitter  because  of 
7.  small  anlge  step  size. 

VIdata=VIdata( : .nzeroindex) ; 

7.  Extract  VQ  data  for  320  frequencies  and  401  azimuthal  angles  at  range  gate  3 

7.  Only  use  2.2  GHz  bandwidth,  corresponding  from  frequency  100  to  frequency  320  (221  freqs). 

7.  VQdata=reshape (data(receiver_pol+l , range_gate ,pol_selection, 100 : 320 , 2 : 402) ,221,401) ; 
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VQdata=reshape(data(receiver_pol+l ,range_gate ,pol_selection, 1 : 320,2 : 402) ,320,401) ; 
7,  VQdata  after  taking  out  repeated  data  due  to  the  pedestal  jitter  because  of 
7,  small  anlge  step  size. 

VQdata=VQdata( : .nzeroindex) ; 


if  pol==’VV’  I  pol==’vv’  I  po^^Vv’  I  pol==  ’  vV1 
I_data=VIdata-ds . avg_of f set_VVI ; 

Q_data=VQdata-ds . avg_of f set_VVQ ; 
elseif  pol==’VHI  I  pol==’vh’  I  pol==’Vh’  I  pol==  ’  vH1 
I_data=VIdata-ds . avg_of f set_VHI ; 

Q_data=VQdata-ds . avg_of f set_VHQ ; 
elseif  pol==’HH’  I  pol==’hh>  I  pol==’Hh’  I  pol==  ’hH’ 
I_data=VIdata-ds . avg_of f set_VHI ; 

Q_data=VQdata-ds . avg_of f set_VHQ ; 
elseif  pol==’HV’  I  pol==’hv>  I  pol==’Hv’  I  pol==  >hV’ 
I_data=VIdata-ds . avg_of f set_VVI ; 

Q_data=VQdata-ds . avg_of f set_VVQ ; 


end 


7.  outdata=sum(I_data. ~2+Q_data. ~2) /221 ;  %  3  April  2003:  Choose  frequencies  from  100  to  320 
7.  outdata=sum(I_data.  ~2+Q_data.  ~2) /101 ;  °i  8  April  2003:  Choose  frequencies  from  190  to  290 
outdata=sum(I_data. "2+Q_data. ~2)/320;  7.  21  April  2003:  Choose  frequencies  from  1  to  320 

ds . theta=theta; 

ds . VIdata=VIdata; 

ds . VQdata=VQdata; 

ds . I_data=I_data; 

ds . Q_data=Q_data; 

ds . IQ_data=I_data+j . *Q_data; 

outdatadB=10*logl0(outdata)  ; 
ds . outdata=outdata; 
ds . outdatadB=outdatadB ; 
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function  Rxhatf ss=estrxhatf s(x, subband) 

7,  Filename:  estrxh.atfs.rn  (Generate  estimated  correlation  matrix  with  forward) 
7,  Author  :  Canh  Ly 
7.  Date  :  21  March  1999 
7. 

7.  Generate  estimated  correlation  matrix  with 
7.  forward  spatial  smoothing  for  scanning  antenna. 

7. 

7.  Syntax:  Rxhatfss=estrxhatfs(x,  subband)  ; 

7. 

7.  where 

/  x:  input  data  matrix 

7.  Rxhatf ss:  estimated  correlation  matrix  with  forward 

7.  spatial  smoothing. 

7. 


mmmmmmmmmfflmramffimmmmxmmm 


7.  This  block  used  each  snapshot  to  compute  the  submatrix  and  averages 
7.  of  all  submatrices. 


7. 


7«P=subband; 

7.  [N,K]=size(x)  ; 

7.M=N-P+1;  7.  Number  of  elements  in  each  subarray 
7.sumsubRx=zeros  (M)  ; 

7«f or  k=l:K 

7.  M=N-P+1;  7.  Reset  new  number  of  length  index  for  each  snapshot 

7.  xtemp=x(  :  ,k)  ; 

7.  for  p=l:P 

7.  xs(  :  ,p)=xtemp(p:M)  ; 

7.  M=M+1 ; 

7.  end 

7.  sum=zeros(size(xs ,  1)  )  ; 

7.  lenxs=size(xs,2) ; 

7.  for  l=l:lenxs 

7.  R=xs(  :  ,l)*xs(  :  ,1)  ’ ;  7.  xs(:,l)’  is  complex  conj  .  transpose 
7.Z=  f  liplr  (flipud(conj  (R)))  ; 

7.  sum=sum+(R); 

7.  end 

7.  tempsum  =  sum/  (P)  ; 

7.  sumsubRx=  sumsubRx+tempsum; 

7.  xs=[]  ;xtemp=[]  ; 

7.end 

7«Rxhatf  ss=sumsubRx/K ; 


nnnrafflfflnmfflranraramfflnnramffl 


mmmmmmmmmmmmmmmiimmmiim 


7.  This  block  uses  the  partition  of  sub-data  matrix  and  then 
7.  computes  the  subcorrelation  matrix  and  then  averages  the 
7.  subcorrelation  matrices. 

7.  K  is  number  of  snapshots,  each  column  of  the  data  matrix,  x. 

7. 
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P=subband ; 

[N,K]=size(x)  ; 

M=N-P+1;  °L  Number  of  elements  in  each,  subarray 

Rxtemp=zeros(M,M) ; 

flag=0; 


for  p=l:P 

Rxtemp=Rxtemp+x(p : M, : )*x(p:M, : ) ’ ; 
M=M+1 ; 

end 

Rxhatf ss=Rxtemp/ (P*K) ; 
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function  [eigval , eigvec] =eigvv(Rxhat) 

7,  Filename:  eigvv.m 
7,  Author  :  Canh  Ly 
7.  Date  :  6/15/95 
7. 

7.  Purpose  :  Computes  eigenvalues  and  eigenvectors 

7. 

7.  Usage  :  [eigval , eigvec]  =eigvv(Rxhat)  ; 

7. 

[V,D]  =eig(Rxhat) ; 
lambdai=diag(D) ; 

[eigval , eigvec]  =  sorteig(lambdai , V) ; 

function  [signal , noise] =snspaces (eigval , eigvec .nsignals) 

7.  Filename:  snspaces.m 
7.  Author  :  Canh  Ly 
7.  Date  :  March  6,  1997 
7. 

7.  Purpose  :  Finds  the  number  of  signal  eigenvalues  and 
7.  theis  indeces  and  splits  noise  and  signal 

7.  space . 

7. 

7.  Usage  :  [signal, noise] =snspaces (eigval , eigvec, nsignals) ; 

7. 

eigvalm  =  eigval(l :nsignals) ; 
eigval index=l : nsignals ; 

7. 

7.  Split  signal  and  noise  space  vectors 

7. 

ncoleigvec=size(eigvec,2) ; 
noiseindex=l ; 
signal index=l ; 

leneigvalindex=length(eigvalindex) ; 
for  j=l :ncoleigvec 

if(j  ==  e igval index (signal index ) ) 
oldsignalindex=signalindex ; 
signal ( : , signal index) =eigvec ( : , j ) ; 
signalindex=signalindex+l ; 
if (signal index>leneigval index) 
signal index=oldsignal index ; 

end 

else 

noise( : ,noiseindex)=eigvec ( : , j) ; 
noiseindex=noiseindex+l ; 
end 

end 
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7,  Filename:  plot  smure  suit  .m 
7,  Author  :  Canh  Ly 

7.  Date  :  25  March  2003 

7. 

7.  Purpose:  Plot  the  SMUSIC  simulation  result. 

7. 

plot (az ,pmusics0, ’b’ ) ; 

7.  set(gca,  ’ Xlim’  ,  [40  58]) 
axis([35  60  0  max(pmusics0)+5] ) ; 
title(plottitlestr) ; 
grid  on 

7.  Filename:  plottgtsetup  ,m 
7.  Author  :  Canh  Ly 

7.  Date  :  13  December  2001 

7. 

7.  Purpose:  Plot  experiment  target  setup. 

7. 

image (A) ; axis  off; hold  off 

7.  Filename:  plotmaveexpdata.m  (Plot  magnitude  average  experimental  data) 
7.  Author  :  Canh  Ly 
7.  Date  :  20  June  2003 
7. 

7.  Purpose:  Plot  magnitude  average  experimental  data 

7. 

plot (ds .theta-1 . 87 ,ds . outdatadB-min(ds . outdatadB) , ’  r ’ ) ; 
axis([35  60  0  max(pmusics0)+5] ) ; 
title(plottitlestr) ; 
grid  on 

7.  Filename:  plotpaveexpdata.m  (Plot  phase  average  experimental  data) 

7.  Author  :  Canh  Ly 
7.  Date  :  20  June  2003 
7. 

7.  Purpose:  Plot  phase  average  experimental  data 

7. 


avecomplex=(l/ size(ds . IQ_data, 1) ) *sum(ds . IQ_data) ; 
plot (ds .theta-1 . 87 , angle (avecomplex) , ’r 1 ) ; 
grid  on 


46 


^function  plotmfreqdata(column) 

7,  Filename:  plotmf reqdata.m 
7,  Author  :  Canh  Ly 
7.  Date  :  20  June  2003 
7. 

7.  Purpose:  Plot  magnitude  data  at  a  particular  frequency. 

7. 

7. 

hold  off 

column  =  get (geo , ’value 1 ) ; 

if  column  ==  size (ds . IQ_data, 1)+1 

plot (ds . theta-1 . 87, ds . outdatadB-min(ds . outdatadB) , ’r ’ ) ; 

else 

plot (ds . theta- 1 . 87, 10*logl0(real (ds . IQ_dat a (column, : ) . ~2+imag(ds . IQ_data (column, : ) . ~2) ) ) ) 

end 

xlabel(’ Azimuth  angle  (deg)’) 
ylabel ( ’ Amplitude ’ ) 


7«function  plotpfreqdata(column) 

7.  Filename:  plotpf reqdata.m 
7.  Author  :  Canh  Ly 
7.  Date  :  20  June  2003 
7. 

7.  Purpose:  Plot  phase  across  scanning  angle  at  a  particular  frequency. 

7. 

7. 

hold  off 

column  =  get (geo, ’value’) ; 

if  column  ==  size (ds . IQ_data, 1)+1 

avecomplex=(l/ size (ds . IQ_data, 1) ) *sum(ds . IQ_data) ; 
plot (ds . theta-1 . 87, angle (avecomplex) , ’r ’ ) ; 

else 

plot (ds . theta- 1 . 87, angle (ds . IQ_dat a (column, : ) ) ) ; 

end 

xlabel(’ Azimuth  angle  (deg)’) 
ylabel ( ’ Phase ’ ) 
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7,  Filename:  plotantpat.m 
7,  Author  :  Canh  Ly 
7.  Date  :  24  March  2003 


7. 

7.  Purpose: 

7. 

7. 

7. 

7. 

7. 

7. 


Plot  5  in  antenna  pattern.  The  data  file  in  c:\smusic_MFRF\mfiles 
was  generated  from  ant=nor5ant;  7.  This  file  output  the  ant  structure 

antraw= [ant . xant  ant . yant] ; 
antsmo= [ant . azscan  ant .magscan] ; 
save  ant5inraw.dat  antraw  -ascii 
save  ant5insmo.dat  antsmo  -ascii 


7. 

7.  25  April  2003 
load  ant46in35GHz . dat 

plot (ant46in35GHz ( : , 1) -3 . 77 , ant46in35GHz( : ,2) , ’k’) 
7.  End  of  25  April  2003 


grid  on 

axis([35  60  0  max(pmusics0)+5] ) ; 
title(plottitlestr) ; 

7.  Shade  the  antenna  beamwidth:  Added  on  21  Mar  2003 
xantf ill= [44 . 7  44.7  51.6  51.6]; 

yantfill=[0  max(pmusics0)+7.5  max(pmusics0)+7.5  0]; 

7.  yantf  ill=  [0  25  25  0]  ; 

tcolor= [0 . 7  0.7  0.7] ; 

p  =  fill (xantf ill, yantf ill.tcolor) ; 

set(p, ’FaceColor’ ,tcolor, ’EdgeColor’ , ’None’ , ’FaceAlpha1 ,0.4) ; 
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^function  plotpfreqant (column) 

7,  Filename:  plotpfreqant  .m 
7,  Author  :  Canh  Ly 
7.  Date  :  21  July  2003 
7. 

7.  Purpose:  Plot  phase  across  scanning  angle  at  a  particular  frequency. 

7. 

hold  off 

column  =  get (geo , ’value 1 ) ; 


if  column  ==  size (ant5 . IQ_data_raw, 2)+l 
load  ant46in35GHz . dat 

plot (ant46in35GHz ( : , 1) -3 . 77, ant46in35GHz( : ,2) , ’k’) 
axis([35  60  0  max(pmusics0)+5] ) ; 

else 


mmmmmmfflmrammmmmmmmiirammrammmmmmmiimmffl 


7.  Old  procedure  of  correcting  phase,  correctantpha  was  created  from  smusic_mfrf _inter_Ver2 .m 

7. 


7.  plot(ant5. theta, unwrap(correctantpha( :  , column))  ,  ’b’) 

7.  hold  on; plot (ant5 . theta, unwrap (angle (ant5 . IQ_data_raw( : , column) ) ) , 1 r 1 ) 
7.  plot  (ant5  .theta,  correctantpha(  :  ,  column)  ,  ’b’ ) 

7.  hold  on;plot (ant5 .theta, angle (ant5 . IQ_data_raw( :, column) ), ’r’ ) 

7. 


mmmmmmmmmmrammmffimmmfflmmmmmxmmmmmm 


7.  New  procedure  of  correcting  phase 

7. 

plot (ant 5 .theta, unwrap (angle (ant 5 . IQ_data_raw_cor ( : , column) ) ) , ’b’ ) 
hold  on; plot (ant5 . theta, unwrap (angle (ant5 . IQ_data_raw_un( : , column) ) ) , ’r ’ ) 


legendstrl= [’ correct  @  freq  ’  num2str (column) ] ; 

legendstr2= [’uncorrect  @  freq  ’  num2str ( column) ] ; 

legend(legendstrl , legendstr2) 

title(’Phase  Profile’) 

grid 

hold  off 

end 


xlabel(’ Azimuth  angle  (deg)’) 
ylabel ( ’Phase  (rad)’) 
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